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Abstract 

A detailed comparison is made between the field-theoretic and ge- 
ometric definitions of topological charge density on the lattice. Their 
renormalizations with respect to continuum are analysed. The defini- 
tion of the topological susceptibility %■> as used in chiral Ward identi- 
ties, is reviewed. After performing the subtractions required by it, the 
different lattice methods yield results in agreement with each other. 
The methods based on cooling and on counting fermionic zero modes 
are also discussed. 
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1 Introduction 



The definition of topological charge density and of topological susceptibility 
on the lattice has by now a long story with contrasting results Jl], |2f. This 
paper intends to be a contribution to clarify the issue. 

Lattice is a regulator of the theory. It should reproduce continuum physics 
in the limit in which the cutoff is removed, i.e., in the limit in which the lattice 
spacing a tends to zero. Like any other regularization scheme, however, 
appropriate renormalizations have to be performed to determine physical 
quantities. Within the rules of renormalization theory the topological charge 
density and its correlation functions can be defined on the lattice with the 
same rigour as for any other operator of the theory. 

In QCD 

Q{x) = ^e, upa F^(x)F; a (x). (1) 

Q(x) has a fundamental physical role, being the anomaly of the Ua{1) singlet 
axial current 

«V» = -2N f Q(x), 

N f _ 

4to = Ev^h^v^). (2) 

1=1 

Nf is the number of light flavours. Eq.(fJ) provides a solution to the Ua(X) 
problem of Gell-mann's quark model in which is conserved and the cor- 
responding £/a(1) is a symmetry, whereas in hadron physics neither parity 
doublets are observed, which would correspond to a Wigner realization, nor 
the inequality < ^/3m n is satisfied, which would correspond to a spon- 
taneous breaking a la Goldstone. 

Eq. (0) could explain the higher value of m v > as suggested by an approach 
based on 1/N C expansion of the theory. At the leading order the anomaly is 
absent being 0(1/N C ), and £7a(1) is a Goldstone symmetry like axial SUa(3). 
The idea behind this expansion is that already at this order the theory de- 
scribes the main physical features of hadrons (e.g. confinement) ||. In the 
1/N C expansion the anomaly acts as a perturbation, displacing the pole of the 
Ua(X) Goldstone boson to the actual mass of the rj '. The prediction is [§, |] 
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where 

X = J d 4 x(0|T(Q(x)Q(0))|0> (4) 

is the topological susceptibility of the vacuum in the unperturbed (iV c = oo) 
theory. This means, among other facts, quenched approximation, fermion 
loops being 0(g 2 N f ) ~ 0(N f /N c ). 

In fact, as we shall discuss in detail below, x in Eq.(f|) is not defined if the 
prescription is not specified for the singularity of the product Q(x)Q(0) as 
x — > 0. In refs. [|J || the prescription which leads to Eq.([|) is the following 

X = J d\x - y)df^(0\T{K ll (x)K v {y))\0) (5) 

where K^x) is the Chern current 

K, = W^S (V* - \9f abc K A l) ( 6 ) 

related to Q(x) by the equation 

d^K^x) = Q(x) (7) 

The prescription Eq.(^|) eliminates all the 5-like singularities in the product 
K )M (x)K u (y) as x — ► y. In any regularization scheme the prescription Eq.(^) 
only leaves a multiplicative renormalization Z 2 for x, Z being the possible 
renormalization of Q(x). Eq. (^) implies that the total topological charge 

Q = J d 4 xQ(x) (8) 

has integer values. 

The regularized version of Q(x), Ql(x), does not obey in general (|7]) [|]. 
According to the general rules of renormalization theory (in pure gauge the- 
ory) 

Q L {x) = ZQ{x). (9) 

In general Z ^ 1 unless Eq. (|7|) is preserved by regularization. To determine 
Z it is sufficient to measure (Ql) = cl 4 Z(Q) on a state belonging to a definite 
eigenvalue of Q (see section 2). 

The product Ql( x )Ql(v) will not satisfy the prescription Eq. (|J) in gen- 
eral at the singularity x — ► y. In the limit a — > it will differ from it 
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by additive terms, 8 X , which can be classified by use of the Wilson operator 
product expansion 0. Defining (V is the 4-volume) 

XL =±Y,QL{x)Q L {y) (10) 

* x V 

we will have 

XL=^Z 2 Q 2 a 4 + 6 X (11) 

where the first term corresponds to the prescription Eq. (|5|). Taking the 
v.e.v. of Eq. (|TTD gives 

XL = a 4 Z 2 x + xo (12) 

with 

Xo=(0\S X \0)- (13) 
Taking the expectation value of Eq. ( |TTD on eigenstates \q n ) of Q gives 

(q n \XL\qn) = ^Z 2 qla 4 + (q n \Sx\q n )- (14) 

It is a generally accepted wisdom that renormalization effects produced by 
short range quantum fluctuations are practically independent of the semi- 
classical instanton background which determines q n . The independence on 
q n of (q n \^x\Qn) can be checked numerically by Eq. ( |T4"D and proves to be 
true within errors ||. Then (q n \8x\<ln) — Xo, and Xo can be determined from 
Eq. ([14] ) as (q n = 0\xL\q n = 0), i.e. as expectation value of xl on the trivial 
topological sector. 
From Eq.([T|) 



X = ^^. (15) 

It is with this prescription that X is expected to be 

X = (180 MeV) 4 (16) 

in the quenched approximation within an 0(1/N C ) systematic error. 

In this paper we will show that, if the prescription Eq.(^) is properly 
implemented, all methods which have been proposed to determine X on the 
lattice give the same result. We shall do this by comparing the geometric 
method | to define Q(x) to the field-theoretic one || 0] for SU(2) 
gauge theory. The same procedure, applied to SU(3) indeed confirms JlTJ 
the expectation Eq. (|T6|) . 
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2 Defining Q(x) on the lattice 



In analogy to any lattice operator, Ql{x) will be defined by the requirement 
that, in the formal (naive) limit a — > 

Q L {x) a ^a 4 Q{x)+0{a & ). (17) 

A prototype definition is 

-1 ±4 

Q L ( x ) = wZ5 51 W CT Tr(n^(x)n pCT (x)). (18) 

tuvpa is the standard Levi-Civita tensor for positive directions while for neg- 
ative ones the relation e^ upa = —i-p UP a holds. 



0(a 6 ) irrelevant terms in Eq . fllTf ) will dissapear in the scaling regime. 
However their presence may be used to improve the operator jT2|. In what 



follows Qi\ (i = 0,1,2) will denote the operator defined by Eq.fll8|) and 
the once and twice improved versions of it respectively. Improvement is 
the recursive smearing of the links developed in |12j| . Also the geometric 



definition Q| com (a;) satisfies Eq. flTTD 



Like any other regularized operator, Ql{x) will mix in the continuum 
limit, when irrelevant terms become unimportant, with all the operators 
having the same quantum numbers and lower or equal dimension. The only 
pseudoscalar of dimension < 4 being Q(x) itself, 

Q L (x) = ZQ(x). (19) 

The naive expectation for Z would be Z = 1 since Q, as an integer, should 
not renormalize. As first realised in |J this is not true on the lattice where 
Ql(x) is not a divergence. Z can be computed in perturbation theory, as it 
was done in the early works on the subject A better way is to measure 
(Ql), Ql being the total topological charge on the lattice, on a state on 
which Q has a known value, e.g. on a one-instanton state where Q = 1. 
This can be done by a heating technique || where a background instanton 
is put by hand on the lattice, and quantum fluctuations at a given value of 
j3 = 2N c /g 2 are added to it. In the continuum the instanton configuration is 
stable, being a minimum of the action, and therefore perturbing it by small 
fluctuations does not change the value of Q. On the lattice instantons are not 
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stable, so that Q could change during this heating procedure. The way to 
avoid this inconvenience is to create a sample of configurations by the usual 
Monte Carlo updating, starting from the original instanton. Each of them is 
checked in its instanton content by a rapid cooling: configurations where the 
original topological charge seems to be changed are discarded. This can be 
done after any number of heating steps and the result must be independent of 
this number. The result is a sample with topological charge Q. Ql measured 
on this ensemble will reach a plateau in this heating procedure, on which Z 
can be read. If the instanton were stable, the plateau would stay flat forever. 

This procedure has been checked and used repeatedly within the field- 
theoretical method [11, [13|, [14], 15[]. The result for Z is shown in Fig. 1 as a 
function of (3 for different definitions of Ql{x). The data for the 0-, 1- and 2- 
smeared operators are taken from ref. |15j and are computed on a 1 -instanton 
configuration on a 16 4 lattice. The data for the geometrical definition are 
new and are computed with the same procedure. For the geometric definition 
Z is compatible with 1 (within two standard deviations). However the values 
of <5| eom have a large spread (as large as Q ± 10) showing that it can assume 
values different from the original Q on configurations which are presumed to 
belong to that sector. Only the average satisfies (Ql) = Q, but not the value 
configuration by configuration. Moreover, on a given configuration the value 
of Q L com depends on the interpolation used to define it 110 . 



3 The topological susceptibility 

The lattice topological susceptibility is written as 

XL = £<Q L (a;)Q L (0)> = ^ (20) 

X * 

and analogously for Q L eom . To make connection to the continuum susceptibil- 
ity as defined by Eq.(^|), in general there will be an additive renormalization 
due to the singularity at x — > y and a multiplicative residual renormaliza- 
tion (x 7^ y) which will simply be the square of Z computed in the previous 
section. 

As a matter of fact (Ql(%)Ql(0)) is expected to be negative due to reflec- 



tion positivity at x ^ since Ql(x) changes its sign under time reversal [|l| 



In fact this holds at distances larger than the extension of the operator if 
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it is smeared. Figs. 2 and 3 show that this is indeed the case both for the 
geometric operator and the field-theoretical definition. Since Q\ is positive, 
its value is determined mainly by the point at x — 0, i.e. by the singularity of 
the product at x — > 0. This peak is there, no matter how Ql(x) is defined, 
and its height depends on the definition used. In Figs. 2 and 3 the values for 
(Ql(x)Ql(O)) have been summed over all points x inside a shell at distance 
\x\ from the origin x — 0. The width of this shell was 1.2 a. 
Thus in general || [7| 

XL = Z((3) 2 a 4 x + M(f3). (21) 

M(j3) will describe a mixing with all scalar operators of dimension < 4 ((3(g) 
is the beta function), 

M((3) = A((3)(^F^F^ + P((3) • 1. (22) 

M is the value of Xo i n Eq.(|T2"D in the lattice regularization. 

To match the prescription of Eq.(|5|), x has to be zero in the sector Q = 0. 
In that sector thus xl = M(j3) and M(f3) can be determined by measuring 
Xl in h. This is again done by a heating procedure ||. The flat, zero field 
configuration, (U^(x) = 1) can be dressed with local quantum fluctuations, 
which do not change its topological content, by the usual updating procedure 
at the desired value of (3. xl will soon reach a plateau: if the sector were 
stable the plateau would persist forever. Instead non-vanishing topological 
charge can be created on the lattice and care must be taken to eliminate 
configurations where this happens. Again this must be done by cooling and 
checks can be done to test the consistency of the procedure. Fig. 4 shows the 
determination of M(f3) for the geometric definition and for three different 
field-theoretical definitions. Analogously to Eq.([T^), from Eq.([2~l~D we obtain 

a\ = XL - m) (23) 
aX Z(0) 2 ' [ ' 

Xl, M and Z depend on the choice of the regulator as well as on the choice of 
the action; a depends on the choice of the action but x must be independent 
of all of it. Fig. 5 shows that this is the case. In this figure we have used 
the data of ref. [H| for the 0- and 2-smeared field-theoretical charges. The 
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data for the geometric definition has been obtained on a 16 4 lattice with 
the same updating procedure (heat bath) and compatible statistics (5000 
configurations). The result of the simulations is in fact xM-i- Usually people 
determine by computing the string tension cr/A| and by assuming the 
physical value for a. This allows to express x 1 ^ i n physical units. We do 
same in order to compare our result with other people's determinations. The 
scale is determined from the data of |17|, . The data at 2-smearings yield 
(x) 1 ^ 4 — 198 ± 2 ± 6 MeV for SU(2) gauge group, the first error being 
statistical and the second one coming from the error in A^. The "naive" 
unsubtracted geometric definition does not scale and is almost one order of 
magnitude larger than the subtracted value. In the jargon of the geometrical 
method this is called an effect of dislocations. It is the mixing to the identity 
operator which indeed describes these dislocations, having dimension lower 
than 4. There is however an additional mixing in Eq. (P2"D which has the same 
dimension as x an d still must be subtracted: checking only by dimension is 
not sufficient to ensure that xl is indeed equal to the physical Xi as defined 
by Eq.®. 

Figs. 6 and 7 show the distribution of values for Ql in the sector with 
trivial topology. Its variance is, apart from a normalization factor, a measure 
of M(f3). A good operator Ql(x) is one for which the subtraction M((3) is 
small compared to xl- On the other hand, also having Z « 1 is more 
reassuring than having a small Z. 

Table I shows xl, Z and M for the 0-, 1-, 2-smeared field-theoretical 
charges and the geometric charge at (3 = 2.57. 

The 2-smeared definition of Ql{x) is the best among these choices. The 
geometric definition is good with respect to Z but is definitively bad with 
respect to the additive renormalization M. 

4 Discussion 

The main conclusion of the above analysis is that with any definition of 
topological charge density on the lattice, an additive renormalization for the 
topological susceptibility and a multiplicative one are necessary. If properly 
renormalized all definitions bring about the same physical value for x- 

Confusion on this subject in the past was generated by a mistreatment 
of renormalization. On the one hand, the geometric definition was believed 
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to be free from renormalizations because it gave always integer values for 
the total topological charge. This seems to be true for the multiplicative 
renormalization. Having integer values, however, does not exempt from hav- 
ing singularities at short distance in the product which defines xl- Fig. 5 is 
clearly proving that. 

The field-theoretical definition started as a naive definition. Z was not 
noticed and put equal to 1; P{0) was subtracted by use of perturbation 
theory. As a result Z 2 x was determined instead of x itself, and found to be 

Hi- 



much smaller than the expectation Eq.(|T( 

The idea was then put forward that the naive definition might not be cor- 
rect and the geometric method || |2(|, the cooling method pT| , |22|] and the 



Atiyah-Singer based methods [23 were developed. The naive method was 
promoted to field-theoretic method only after introducing Z and a correct 
subtraction M [BL 0]. The non-perturbative determinations of these con- 
stants P] as explained above, finally brought about a reliable determination 
of Xi which is indeed regulator independent. 

The cooling method authomatically performs the additive subtraction be- 
cause it gives xl = on the trivial sector; and also brings Z to 1 by freezing 
the quantum fluctuations. The problem with it was that instantons could 
be lost in the procedure, leading to an underestimation of x- Cooling with 
improved forms of the action |23, |25| seems to have eliminated this problem 
and gives indeed results which confirm the field-theoretic determination. The 
same seems to be true for the modern versions of the Atiyah-Singer proce- 
dure |26| . 
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Figure captions 



Figure 1. Values of Z as obtained by the heating method for the geometric and 0- 
, 1- and 2-smeared field-theoretic topological charges, (circles, squares, 
up-triangles and down-triangles respectively). 

Figure 2. Correlation function (Ql(x)Q l (0)) as a function of |x|/(1.2 a) for the 
0-, 1- and 2-smeared topological charges (squares, up-triangles and 
down-triangles respectively) at f3 — 2.57. 

Figure 3. The same as in Fig. 2 for the geometrical topological charge. 

Figure 4. Values of M as obtained by the heating method for the geometric and 0- 
, 1- and 2-smeared field-theoretic topological charges, (circles, squares, 
up-triangles and down-triangles respectively). 

Figure 5. x i n A-l and MeV units for the unsubtracted geometrical charge (stars), 
subtracted geometrical charge (circles) and 0- and 2-smeared charge 
(up and down triangles). 

Figure 6. Distribution of Ql in the zero-topological charge sector Q = for the 
0-smeared (solid line) and 2-smeared (dotted line) topological charge 
densities at f3 — 2.57. 

Figure 7. The same as in Fig. 6 for the geometrical topological charge. 



Table caption 

Table I. xl, Z and M for the 0-, 1-, 2-smeared and geometric topological 
charge density operators at (5 = 2.57. 
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Table I 



operator 


10 5 x XL 


Z 


10 5 x M 


geometric 


16.6(3) 


0.937(26) 


13.26(23) 


0-smeared 


2.320(52) 


0.240(26) 


2.200(32) 


1-smeared 


1.010(49) 


0.507(9) 


0.440(18) 


2-smeared 


1.165(64) 


0.675(8) 


0.187(5) 
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Fig. 3 
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